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Abstract 



We study effects of local electron interactions on the persistent current 

of one dimensional disordered rings. For different realizations of disorder 

we compute the current as a function of Aharonov-Bohm flux to zeroth and 

first orders in the Hubbard interaction. We find that the persistent current is 

enhanced by onsite interactions. Using an avoided level crossings approach, we 

derive analytic formulas which explain the numerical results at weak disorder. 

The same approach also explains the opposite effect (suppression) found for 

spinless fermion models with intersite interactions. 
72.10.-d, 71.27.+a, 72.15.Rn 
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I. INTRODUCTION 



Recent experiments have found that small isolated metallic rings, threaded by magnetic 
flux, carry persistent currents. Although the existence of this effect had been anticipated 
theoretically for many years Jlj-Q, the magnitude of the observed current was found to be 
much larger than expected. 

Experiments performed on moderately disordered rings with small transverse width have 
measured the magnetic response of an ensemble of copper rings 0, of isolated gold loops 
0, and of a clean, almost one-dimensional, sample [[7|. 

The ensemble averaged persistent current was found to have an amplitude of order 
10~^et'//L and periodicity $o/2, (where L is the circumference of the ring, vj is the Fermi 
velocity, e is the electron charge, and $o = h/e is the flux quantum). The single loops 
experiment reported flnding an unexpectedly large current of order eVf/L, the free electron 
value, and displaying periodicity in flux with period $0- 

Theoretical approaches can be classifled into non-interacting and interacting electrons 
approaches. Non interacting theories which take into account the effects of disorder and 
deviations from a perfect ID geometry predict values of the average persistent current 
that are smaller than observed, but still within one or two orders of magnitude of the ex- 
perimentally determined value of the average current . More surprisingly, the theoretical 
estimates for the typical value of the persistent current, measured in the single loop exper- 
iment, are perhaps two or three orders of magnitude smaller than the reported values. It 
seems that agreement between theory and experiment worsens with increasing disorder. This 
suggests that perhaps interaction corrections are quantitatively important in the moderately 
disordered regime. 

A flrst guess would be that interactions enhance the impurity scattering and further 
suppress the persistent current, as happens for the conductance of the disordered luttinger 
liquid . However, this intuitive analogy may be quite misleading. While the conductance 
depends on the Fermi surface properties, i.e. velocity and scattering rates, the persistent 
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current is a sum of contributions from all occupied states. Also, by Galilean invariance, an 
interacting electron liquid without disorder carries the same current as the non interacting 
gas, even though its Fermi velocity may be renormalized. Thus it is plausible that in the 
disordered system, interactions could work both ways, i.e. either suppress or enhance the 
persistent currents under conditions which need to be explored. 

Exact diagonalization studies of spinless fermions with intersite interactions on small 
rings [0,0 , have found that weak interactions reduce the persistent current below its non- 



interacting value. On the other hand, recent reports [0,0 for Hubbard interactions found 
the opposite effect. 



In this paper [T£| we clarify this seeming contradiction, and gain some insight into the 
puzzle of interaction corrections. In Section |I| we diagonalize the tight binding Hamiltonian 
numerically for different realizations of disorder, and compute the persistent current to 
zeroth and to first order in the Hubbard interactions. We find that interactions enhance 
the persistent current. In Section |T| we present the Avoided Levels Crossing (ALC) theory, 
which provides an analytic approximation to the numerical results. This is an expansion, 
at weak disorder, of nearly degenerate eigenstates at fluxes close to points of time reversal 
symmetry. The correlation between the zeroth and first order currents is explained by a 
common mechanism, i.e. the avoided level crossings. The opposite (suppression) effect 
for interacting spinless fermions, is also explained by the ALC theory in Section We 
conclude by a summary and future directions. 

II. NUMERICAL PERTURBATION THEORY 

We consider a tight-binding Hamiltonian on a periodic chain with repulsive on-site Hub- 
bard interaction: 

H = Ho + UY^n^^m^, (1) 

i 

Ho = -^[e*^^^^a|+i^ai^ + h.c.] + ^e^atoi^. (2) 



where a}^ creates an electron at site i with spin a. Our unit of energy is the hopping 
energy, is the dimensionless on-site disorder energy uniformly distributed in the interval 
[—W/2,W/2]. The chain has L sites, Ng electrons, the flux through its center is given by 
With these conventions, the energy spectrum is periodic in the enclosed flux with period 
$0) s-nd the current is given by 

m - --^m, (3) 

where 

E{^) = Eoi^) + Ei{<5) + 0{U^). (4) 

Eq is the exact ground state energy of Hq, whose single electron eigenstates are determined 
numerically. Ei is the first order correction in U, which is given by 

El = Uj2ni^nii 

i 

riia = (^olalcraicrl^o) (5) 

where ^'o is the Fock ground state of Hq. Thus by diagonalizing Hq we can readily obtain 

/($) = /o($) + /i($) + 0(^/2). (6) 

The effects of disorder on /o($) can be seen in Fig. 1, where it is shown as a function of 
flux for a half filled lattice of six sites, for different disorder strengths W. Disorder smoothens 
and reduces the magnitude of Iq, and for >> 1 it is dominated by the first harmonic 

sin(27r$/$o). 

In Fig. 2, is plotted. There are several features of the first order interaction 

correction that deserve comment. First, the most important observation is that it generally 
enhances the noninter acting current for all values of disorder strength and flux. That is to 
say: there is a positive correlation between the non interacting current and the first order 
correction, 

(/o($)/i($)),,$ >0. (7) 
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This result is found numerically to hold for all realizations of disorder which we have used. 
Second, in the limit of weak disorder, becomes singular at discontinuity points of 

/o($), which are at fluxes (m + |)$0; ('^'^'o) for even (odd) number of filled orbitals, where 
m are integers. Finally, we see the first order current also becomes dominated by its first 
harmonic at large disorder. 

To study the scaling properties of the current with system size and disorder strength, 
it is convenient to characterize the strength of both Jq and Ji by their amplitude at $ = 
$o/4. And while it does not capture the singular nature of near the regime of very 

weak disorder, this characterization is still useful to study the scaling properties of the 
first order current, even in this limit. Previous studies have found [§,0] both numerically 
and analytically, that the amplitude of the noninteracting current, averaged over disorder, 
behaves like: 

|/o($o/4)|^^^exp(-L/0. (8) 

By fitting to this functional form we extracted the behavior of the localization length ^ 
for different system sizes and disorder values. Figure |^ shows the values of the localization 
length obtained for sizes L=6,10,14,20 and 25, at half-filling, averaged over many (up to four 
thousand) impurity configurations for each value of W. We find good agreement between 
our inferred value of the localization length and the known asymptotic form in the weak 
disorder limit: ^ = 105/W'^, valid when W 27t for large systems at half-filling. 

There is less agreement in the strongly localized limit where the localization length is 
known to behave like: ^ = (ln(W^/2)-l)-^ We find a better fit taking ^ = 1/ \n{W°-^ /2.5) for 
the L = 25 data. This discrepancy could be due to the small sizes considered, or a breakdown 
of Eq. (H) when ^ is of order unity. We have calculated the amplitude /i($ = $o/4) as 
function of strength of disorder, characterized by the scaling parameter L/^, for system sizes 
of L=6,10,14,20,25. At weak disorder (large ^), the amplitude increases with the strength 
of disorder achieving its maximum value at some intermediate strength of disorder ( ^ ~ L 
). In this weakly disordered regime, a single scaling function could be used to describe the 
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results: 

|Ji|~/(L/0/L. (9) 

This behavior is similar to that of the amplitude of the noninteracting current given by (^. 
Upon increasing the impurity scattering further, the first order current decreases with dis- 
order. When the single particle wavefunctions are sufficiently localized, Ii can be described 
by a different scaling form: 

\Ii\r^giL/0, (10) 

where g is a. decreasing function of its argument. Eq. (|T0|) suggests that in the localized 
regime Ji dominates Jq for large system sizes. However for localized doubly occupied single 
electron states 

{i^a]i^al\Uj2^^^^il\ (11) 

Thus even for weak interactions, the interaction corrections may be larger than the non 
interacting level spacings which go as 1/L. This invalidates perturbation theory in U in the 
localized regime. 

III. AVOIDED LEVEL CROSSINGS THEORY 

Here we will discuss how the numerical results of the previous section can be understood 
in terms of avoided level crossings (ALC) at weak disorder. First, we derive the ALC 
approximation for the noninteracting current Jq. 

A. ALC theory for Iq 

In Fig. I we can see a typical spectrum for a tight-binding Hamiltonian of a 6 site ring, 
as a function of the applied flux. In the absence of disorder {W = 0) the eigenenergies are 

6„(<l>) = -2cos(^(n+|-)), (12) 
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with period $o = 27r. At points of time reversal symmetry, i.e. where $ is an integer 
multiple of $o/2, level crossings occurs between states of opposite angular momenta. The 
noninteracting persistent current is a sum over the currents carried by all occupied levels, 

a p occupied o 

g 1,^*). (13) 

where n,s are the orbital and spin index, respectively. The noninteracting current /o($) will 
be a smooth function of $ away from the points of level crossings. By symmetry of $ — — 
any pair of levels cross with opposite slopes, and thus if they are both fully occupied their 
contribution to the total current cancels. The only nonvanishing contribution comes from 
a topmost level which is not compensated by its partner. Then the current changes sign 
abrubptly as the occupation moves from one branch to another. Thus for an odd number 
of fully occupied orbitals the current discontinuity occur at $ = (m + l/2)$0; otherwise 
the discontinuities will occur at $ = m$o- In between the discontinuities the current varies 
linearly with the flux, which explains the periodic sawtooth shape for Jq as seen in Fig. |1|. 
Introducing a small amount of disorder lifts the degeneracy at the crossings by opening small 
gaps. This reduces the persistent current Jq since the occupied levels have smaller slopes 
near the former crossing points. This explains the behavior shown in Fig. ^ weak impurity 
scattering softens the discontinuities and leads to an overall reduction of the magnitude of 
current. We can quantify this observation by examining pairs of levels with momenta ±k 
which cross at $ = 0. We specialize to the case of an even number of filled levels with equal 
occupations for both spin directions. The unperturbed energies are given by 

/ 27r $ \ 

e±fc($) = -2cos ±/. + -— (14) 



L $0. 

with k = 271771/ L, for a positive integer m. Consider the effect of a weak random potential 
ej which lifts the degeneracy in the 2x2 subspace of k, —k, 



Ha 



(15) 



where 
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k,-k 



^exp(-i2fca;j)ej. 



1=1 



From now on we will omit the subcripts /c, —k off V . 



The eigenvalues and normalized eigenfunctions of ([15D are: 



where: 



A. 



B. 



1 + 



4IVI 



efc-e-fc + \/(e-fc-efc)2 + 4|V|2 
-2\/* 



-1/2 



and. 



which satisfy: 



6+($) + e_($) = efc($)+e_fc($). 



(16) 



(17) 



(18) 



(19) 



(20) 



We see that the contribution of the occupied orbitals to Jq is unchanged by weak disorder 
since: 



jk, — k 
'O 



(21) 



Let us now use (|T3]) to calculate /o(^') for free electrons with weak disorder by summing 
over occupied levels. We fill all levels up to A;/ = 27rm//L, plus two electrons at kj so 
that an even number of orbital levels is filled. The total noninteracting current is given by 
the contribution from the filled level pairs, together with the contribution of lowest m = 
(nondegenerate) level and the topmost orbital level at m/: 
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mt — l 



-2 E 



47r 2tt / $ 
— — sm — m + — 
L$n L \ $0 



-(mf—1) u 

Neglecting corrections of order we obtain the expression: 

Svr sin(A;/) 



9$ 



(22) 



/o($) 



$0 V$o 



L\V\ 



47r sin(A;/) 



(23) 



for $ G [— $o/2, $o/2]- This expression for Iq is vahd as long as the energy scale of the 
disorder is much smaller than the level spacing at the Fermi level: 



\V\ « A'Ksin{kf)/L. 



(24) 



Using the relation between \V\ and the mean free path /ez defined by the one dimensional 
Born approximation, 



271 sin^ k 



el 



f 



(25) 



According to (0), the ALC approximation is valid for 

8rtlei/L > 1, 

which is the ballistic, or delocalized regime.^] In terms of l^i one can write 

87rsin(fc/) 



(26) 



/o($) 



I 



$0 



+ 



L 

8^ 



el 



(27) 



In order to compare ( pTf ) to the numerical result for Eq. (0), one needs to determine the 
parameter Igi. Since fluctuations in Jq for different disorder realizations, we determine lei by 
fitting (^) to the numerical Iq. Once lei is determined for a particular disorder realization, 
we use it to evaluate Ji as shown below. 



^There is no intermediate "diffusive regime" between the localized and ballistic regimes in one 
dimension 
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B. ALC theory for h 

We shall now proceed to use the same approximation to explain the behavior of 
According to (||), the first order energy is a sum of the density squared at all sites, 

E^{^) = UY.n,^n,^ (28) 

i 

where for the filling of an even number of orbitals per spin, the single spin density is, 

mj—l 

n,,($) = |^^o(0r + l^™„-(Or + E E (29) 

m=l (T=it 

We will now show that the first order energy, 

E,{^) = UY.n,^n,^ (30) 

i 

is enhanced near level crossings. 

For a pair of fully occupied levels one has, by unitarity 

+ \i^m,-{i)? = \M^)\' + \i^-ki^)\' = \. (31) 

Consequently, for weak disorder, all occupied levels, except for the last one, contribute 
constant values to the total density: 

n.T = ^ + l^-«r (32) 
or in terms of the coefficients A, B (O) 



The first order energy is thus given by (|32D , 



EmiU=^^^+'^^^. (34) 



^This approach is similar to the discussion of Aharonov-Bohm oscillations of the participation 



ratio in disordered rings Ref. |16] 
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Using ([T^) and (pSj) we can write 



^ (2mf + 1)2 1 



1 + 



21 -1 



(35) 



Differentiating (|35|) with respect to flux yields 



Ji($)/f/ 



(L$o)2 



$0 



212 



(36) 



In Fig. 1^, (|3^) is compared to the numerical result for Ii for various values of disorder. 
For each disorder realization, l^i is determined by fltting the numerical and ALC results 
for Jq. We see that for weak disorder there is a satisfying agreement between the ALC 
approximation and the numerical results. In Fig. the disorder is too large, and the ALC 
approximation fairs badly. 

Eq. (^) explains both the positive correlation between Iq and Ii of (|^). Since in one 
dimension the mean free path (lei) and the localization length (^) differ by a proportionality 
factor of order unity ||l5l, (|36D agrees with the empirical scaling form (||) ~ f{L/lei)/L 
which was found numerically at weak disorder. 



IV. DIFFERENCE OF HUBBARD MODEL AND SPINLESS FERMIONS 



In recent papers |TO|,|TT|, small disordered rings of spinless electrons have been exactly 
diagonalized. The persistent currents have been computed as a function of disorder and 
interaction strength. In contrast to our result for the Hubbard model, interactions have were 
seen to reduce the persistent current at weak disorder. Here we apply the ALC approach to 
explain this apparent difference between the models. 

The spinless Hamiltonian H"^ is given by 

i 

^J = -E[e'^^^^4+i«. + h.c.] + ^Qafa,. (37) 

i i 

The non interacting current, Jg, is given by half the value of Eq. ([TB|). In the ALC approx- 
imation, it is given by half the value of Eq. (|27|). 
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Following the analogous derivation of Ji, we use (^) to obtain: 



where 



E'^{<l>) = (2m/ + 1)^^(0) + 2\A^\^\B^\^U{2kf), 



which implies 



El{<!>) = (2m/ + l)2f/(0) + U{2kf) 



1 + 



SttL 



L J V$o 

$0 



For the spinless nearest neighbor case studied in Ref. [ITT 



Ei($) = UY^Uim+i, 

i 

the corresponding Fourier coefficient is 

Ui2kf) = 2Ucos{2kf). 



(38) 



(39) 



(40) 



(41) 



(42) 



Above a quarter filling, kf > 7r/4, U is negative. Consequently I'l has the opposite sign 
to that of Ji in Eq. (0), and to Jg. That is to say, the persistent current of the spinless 
fermion model is suppressed by the interactions. 

The difference between the models can be attributed to the different effects of intersite 
interactions (in the spinless model), versus local interactions (in the Hubbard model). 



V. REMARKS AND CONCLUSIONS 

We have investigated the first order effect of Hubbard interactions on the persistent 
current in one dimensional disordered rings. 

The findings can be summarized as follows. Ji was found to correlate in sign with Jq, a 
fact which goes contrary to the naive intuition that interactions suppress currents (e.g. as 
they do for the conductivity of the Luttinger model with an impurity P|). 
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We can understand this result by observing that /i depends on charge fluctuations which 
vary strongly with flux near the degeneracy points, which also determine the sign and 
magnitude of Jq. Using the avoided level crossings theory, we obtain analytical expressions 
which flt the numerical results for both Jq and Ji. We can use the same theory to explain 
why spinless fermions with intersite interactions exhibit suppression of currents rather than 
enhancement. 

We have found numerically that Ii{L/^) scales differently with the localization length 
than Iq{L/^), it seems that effects of electron-electron interactions grow as disorder is in- 
creased. 

However, first order perturbation theory in the weakly disordered case is expected to hold 
as long as the interaction matrix elements do not exceed the single particle level spacings. 
Thus we are restricted to the regime U < 1. We can draw on the exact diagonalization 



results |1T0|JTT[] where for weak disorder, the numerical currents are found to vary linearly 
with interaction strength in a sizeable regime. Thus we believe that first order perturbation 
theory should be valid for physically interesting interaction parameters. 

It would be very satisfying if a similar analysis could be applied to the experimentally rel- 
evant case of three dimensional rings. Preliminary results yield positive correlations between 
Jo and Ji, but where a simple minded application of the ALC approach cannot describe the 
diffusive regime of l^i « L. It would also be important to understand the effects of true 
long range Coulomb interactions with the screening and exchange effects which are absent 
in the Hubbard model. 
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FIGURES 

FIG. 1. The noninteracting current (Iq) as a function of the applied flux for several strengths 
of disordered potential: IF=0, 0.5, 1.0, 2.0 and 4.0. All curves are for a half-filled lattice of six 
sites. 

FIG. 2. The first order interacting current (/i) as a function of the applied flux for disordered 
potential values: W= 1.0, 2.0 and 4.0. As before, all curves are for a half-filled lattice of six sites. 

FIG. 3. The logarithm of the localization length (In^), implied by Eq. (|8|), is plotted as a 
function of strength of disorder (W) for half-filled rings with A^=6, 10, 14, 20 and 25 sites. Data 
for the different system sizes collapse well onto a single curve. 



FIG. 4. Energy level spectrum of the tight-binding IIamiltonian(12) as a function of applied 



flux for a A^=6 ring in the absence of disorder (solid lines) and for weak disorder (dotted lines). 

FIG. 5. Comparison between numerically determined interactions correction Ii (solid lines), 
and the analytic ALC result ( [3^ ) (dashed lines). Figs, (a) - (c) show results for three values of 
disorder strength W. lei{W) are determined by fitting Eq. ( p7| ) to the numerical disorder averaged 
/o($). Error bars depict fluctuations of numerical Ii for different disorder realizations. 
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